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ABSTRACT 

The  viscous  flow  region  in  low  density  hypersonic  ax i symmetric 
nozzles  was  investigated  both  theoretically  and  experimentally.  Non- 
similar  solutions  were  obtained  for  the  internal  laminar  boundary  layer 
equations  which  include  second-order  transverse  curvature  terms.  These 
solutions  were  obtained  on  a  CDC  1604  digital  computer.  Four  existing 
low  density  axisymmetric  nozzles  were  considered.  Numerical  solutions 
were  obtained  for  these  nozzles  using  various  plenum  chamber  conditions 
and  wall  temperature  distributions.  The  plenum  chamber  conditions  used 
in  the  numerical  solutions  for  these  nozzles  produced  a  test  section 
Mach  number  range  of  about  3.0  to  18.0  and  a  test  section  Reynolds 
number  per  foot  range  of  about  1000  to  15,000.  Some  results  of  the 
numerical  solutions  are  compared  with  experimental  measurements  of 
pitot  pressure,  relative  heat  flux,  and  nozzle  wall  heat  transfer 
coefficients.  In  all  comparisons  the  numerical  solutions  appear  to  be 
consistent  with  the  experimental  data.  The  results  of  this  investi¬ 
gation  indicated  that  the  boundary  layer  equations  adequately  described 
the  viscous  region  in  a  nozzle  where  the  mean  free  path  Is  of  the  order 
of  one-tenth  of  an  inch;  and  furthermore,  these  equations  adequately 
described  the  viscous  region  where  the  boundary  layer  thickness  was  of 
the  order  of  the  nozzle  radius. 
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NOMENCLATURE 

py/peMe 

Constant  pressure  specific  heat 

Skin  friction  coefficient 

Defined  in  Equation  22 

Enthalpy  ratio,  H/Hc 

Local  total  specific  enthalpy 

Total  specific  enthalpy  in  plenum  chamber 

Static  specific  enthalpy 

Thermal  conductivity 

Reference  body  length 

Mach  number 

Total  pressure 

Static  pressure 

Pitot  pressure 

Prandtl  number,  uCp/k 

Nozzle  wall  heat  transfer 

Radius  of  curvature  of  converging  portion  of  nozzle 

Longitudinal  radius  of  curvature  in  meridian  plane 

Reynolds  number,  pux/y 

Nozzle  wall  exit  radius 

Nozzle  wal I  radius 

Defined  by  Equation  4 
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Static  temperature 
Total  temperature 

Total  temperature  behind  a  normal  shock 
Adiabatic  wall  temperature 

Transverse  curvature  term  defined  by  Equation  31 

Velocity  component  in  x  direction 

Velocity  component  In  y  direction 

Distance  along  nozzle  wall 

Distance  norma  I  to  nozzle  wall 

Total  axial  distance 

Distance  along  nozzle  axis 

Nozzle  wal I  angle 

Dimensionless  velocity  gradient,  see  Equation  27 
Boundary  layer  thickness  defined  as  value  of  y  where 
u/ue  =  0.995 
Displacement  thickness 
y/« 

Transformed  y  coordinate 
Momentum  thickness 

Heat  transfer  coefficient  defined  by  Equation  55 
Viscosity 

Transformed  x  coordinate 
Mass  density 


Stream  function 
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Subscripts 

(£  Nozzle  centerline 

e  Edge  of  isentropic  core 

w  Nozzle  wal I 
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SECTION  I 
INTRODUCTION 

Since  the  advent  of  flying  vehicles  there  has  existed  the  need 
to  simulate  for  aerodynamic  testing  purposes  the  environment  in  which 
these  vehicles  travel.  With  the  recent  reality  of  vehicles  traveling 
in  the  upper  atmosphere  of  earth,  and  possibly  other  planets,  this  need 
has  become  more  important  due  to  the  complication  of  problems  associ¬ 
ated  with  an  unfamiliar  environment. 

The  flow  field  about  a  vehicle  moving  In  the  earth’s  atmosphere 
at  altitudes  of  about  200,000  feet  to  400,000  feet,  is  characterized  by 
low  Reynolds  number  and  usually  high  Mach  number.  As  is  well  known,  when 
the  Reynolds  number  is  low  the  boundary  layer  is  laminar  and  thick. 
Therefore,  an  inherent  problem  In  the  design  of  a  nozzle  for  a  wind 
tunnel  which  is  to  simulate  this  flow  regime  is  the  analysis  of  thick 
laminar  boundary  layers.  This  can  be  exemplified  by  considering  an 
existing  low  density  Mach  ten  nozzle,  which  will  be  discussed  later, 
that  has  been  used  to  take  much  useful  aerodynamic  data.  The  boundary 
layer  thickness  Is  about  80  per  cent  of  the  nozzle  radius  which  means 
that  96  per  cent  of  the  actual  area  of  the  nozzle  is  dominated  by 
viscous' effects.  The  necessity  of  Investigating  boundary  layers  in 
low  density  hypersonic  nozzles  is,  therefore,  rather  obvious. 

The  purpose  of  this  investigation  is  to  theoretically  study 
the  viscous  flow  field  associated  with  low  density  hypersonic 
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ax i symmetric  nozzles,  and  to  compare  experimental  data  with  the  results 
of  the  theoretical  solutions. 

The  approach  taken  here  is  to  obtain  general  non-similar  solu¬ 
tions  to  an  appropriate  set  of  boundary  layer  equations  which  contain 
second-order  transverse  curvature  terms.  These  equations  are  trans¬ 
formed  into  the  Sn-plane  for  convenience  of  numerical  solution.  A 
general  discussion  of  the  procedure  used  to  obtain  the  numerical  solu¬ 
tions  is  given.  Solutions  for  various  conditions  in  four  low  density 
nozzles  will  be  presented,  followed  by  a  comparison  of  theoretical 
results  and  experimental  data. 

Three  of  the  four  nozzles  considered  use  arc  hooters  to 
generate  high  enthalpy  conditions  in  the  plenum  chamber.  These  high 
enthalpy  conditions  force  one  to  investigate  the  possibility  of 
ionization  and  dissociation  due  to  the  high  temperatures  produced  by 
the  arc.  Since  nitrogen  is  the  test  gas  used  in  each  nozzle,  it  is 
possible  to  use  a  plenum  chamber  which  is  sufficiently  large  tn  permit 
recombination  and  also  sufficiently  small  to  prevent  an  unreasonable 
amount  of  energy  loss  due  to  heat  transfer.  The  analysis  is  therefore 
simplified  since  the  specie  concentration  equations  do  not  have  to  be 
included  In  the  governing  system  of  equations. 

The  design  of  low  density  hypersonic  nozzles  is  usually 
accomplished  by  applying  a  boundary  layer  correction  to  some  desired 
inviscid  flow.  One  such  method  is  due  to  Potter  and  Durand  (I)  ,  which 
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was  recently  modified  by  Potter  and  Carden  (2).  This  method  is  based  on 
an  integral  technique  which  uses  the  similar  solutions  of  Cohen  and 
Reshot ko  (3).  The  results  given  by  this  method  have  proved  to  be  very 
useful  in  designing  a  nozzle  for  a  particular  test  condition.  On  the 
subject  of  non-similar  Internal  laminar  boundary  layers  the  literature 
appears  to  be  rather  limited.  To  the  author's  knowledge  there  presently 
exist  no  non-similar  solutions  to  the  internal  laminar  boundary  layer 
equations  for  a  specified  nozzle  geometry.  This  situation  is  somewhat 
unfortunate  since  it  is  desirable  from  both  the  practical  and  economic 
point  of  view  to  be  able  to  determine  what  flow  conditions  might  be 
expected  from  one  nozzle.  Similar  solutions  are  not  applicable  for  such 
an  investigation  since  similarity  cannot  be  satisfied  in  general  for 
specified  plenum  chamber  conditions,  nozzle  geometry,  and  wall  temper¬ 
ature  distribution  all  occurring  simultaneously.  The  present  work  was 
undertaken  in  order  to  provide  a  means  of  analyzing  nozzles  without  the 
restriction  of  simi lar  boundary  layers. 
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SECTION  II 

MATHEMATICAL  ANALYSIS 

In  order  to  take  maximum  advantage  of  the  integration  tech¬ 
niques  developed  by  Jaffe,  Lind,  and  Smith  (4),  the  mathematical  for¬ 
mulation  of  this  internal  flow  problem  is  handled  similarly  to  their 
work  for  external  flow.  The  governing  system  of  equations  is  trans¬ 
formed  from  the  physical  xy-plane  into  a  Cn-plane  by  appropriate 
transformation  of  variables.  Certain  boundary  layer  parameters  are 
derived  and  the  numerical  solution  of  the  transformed  equations  is  dis¬ 
cussed. 

2.1  GOVERNING  SYSTEM  OF  EQUATIONS 

The  governing  equations  are  taken  from  Probstein  and  Elliott 
(5).  These  equations  In  curvilinear  coordinates  including  second-order 
transverse  curvature  terms 'are: 

continuity  equation 


3(pru)  3(prv)  _  n 
3x  3y 


momentum  equation 


pu 


pv 


a  u  _  dp 

8y  dx* 


a 

ay 


( I ) 


(2) 
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static  energy  equation 


3h 

pu  ^  +  pv 


«*.  =  udE  +  ii- 

3y  u  dx  r  3y 


V  on  ou  l 

r  PF  3?  +  v  s?J 


(3) 


The  coordinate  system  is  defined  In  Figure  I  with  the  r(x,y)  term  being 
defined  for  internal  flow  as 


r(x,y)  =  rQ(x)  -  y  cos  a  (4) 

Probstein  and  Elliott  obtained  Equations  I,  2,  and  3  by  the 
usual  order  of  magnitude  analysis  of  the  general  form  of  the  continuity 
equation,  Navler-S+okes  equations,  and  energy  equation  expressed  in 
curvilinear  coordinates.  The  assumptions  made  in  the  analysis  were  that 
6/Rl  Is  small  compared  to  unity  and  6/r0  is  of  the  order  of  unity. 

Since  the  boundary  layer  thickness  may  be  of  the  order  of  the  nozzle 
radius.  Equations  I,  2,  and  3  are  valid  for  nozzles  which  have  a  longi¬ 
tudinal  radius  of  curvature,  R|_,  much  larger  than  the  nozzle  radius,  rQ. 
This  stipulation  is  normally  satisfied  in  axlsymmetrlc  nozzles  which 
are  used  in  low  density  wind  tunnels. 

The  axl symmetric  boundary  layer  equations,  which  contain 
second-order  transverse  curvature  terms,  can  be  obtained  from  Equa¬ 
tions  I,  2,  and  3  by  replacing  r  with  rQ.  Since  rQ  Is  a  function  of  x 
only,  it  is  eliminated  from  Equations  2  and  3  and  therefore  only  appears 
in  the  continuity  equation.  The' resulting  axisymmetrlc  equations  can 
be  used  to  describe  internal  or  external  boundary  layers. 

Equations  I,  2,  and  3  differ  slightly  from  the  corresponding 
equations  of  Probstein  and  Elliott  In  that  the  specific  heat  and  Prandtl 
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number  are  not  necessarily  assumed  to  be  constants.  However,  all  numer¬ 
ical  results  presented  herein  are  for  constant  Prandtl  number  and  con¬ 
stant  specific  heat. 

A  more  convenient  form  of  the  energy  equation  for  this  work  is 
the  total  energy  equation.  It  is  derived  by  multiplying  the  momentum 
equation  by  u  and  adding  the  result  to  the  static  energy  equation  to 
give 


3H 

pU  37 


3H 
Pv  3y 


i  !_ 

r  3y 


r 


L  L 


vi 

Pr 


3H 

3y 


~  +  V 


I 

F r 


where 


H  = 


h  + 


u£ 

2 


(5) 


(6) 


As  noted  by  the  total  derivative  of  p  with  respect  to  x  in 
Equations  2  and  3,  the  y  component  of  the  momentum  equation  is  given  by 


- 1£.  -  0  (7) 

5y 

The  validity  of  this  equation  is  sometimes  questioned  for  very  thick 
laminar  boundary  layers.  However,  the  analysis  of  Probstein  and  Elliott 
indicates  that  this  equation  is  consistent  with  the  other  equations  in 
the  governing  set.  Also,  the  numerical  solutions  of  the  governing 
equations,  which  implicitly  contain  Equation  7,  are  consistent  with 
experimental  data. 

Expressing  yas  an  arbitrary  function  of  T 

v  =  vCT)  (8) 
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and  taking  the  thermal  equation  of  state  as 

p  =  pRT  (9) 

one  has,  with  Equations  I,  2,  and  5  and  the  appropriate  boundary  con¬ 
ditions,  the  governing  system  of  equations  whose  solution  is  desired. 

For  the  boundary  conditions  at  the  edge  of  the  boundary  layer, 
it  is  assumed  that  an  Isentropic  core  flow  exists  along  the  centerline 
of  the  nozzle.  There  is  no  restriction  on  the  size  of  this  core,  as 
long  as  the  gas  properties  along  the  centerline  can  be  found  from  the 
isentropic  flow  relations.  The  boundary  conditions  at  the  nozzle  wall 
are  taken  as  zero  velocity  and  a  prescribed  wall  temperature  distribu¬ 
tion.  A  specified  wall  heat  transfer  distribution  could  be  used  in 
place  of  the  wall  temperature  distribution;  however,  the  boundary  con¬ 
dition  used  herein  is  the  latter. 

2.2  TRANSFORMATION  OF  THE  GOVERNING  EQUATIONS 

2.2.1  Transformation  Variables 

The  transformation  variables  used  are  a  combination  of  a 
modified  form  of  the  Mangier  transformation  proposed  by  Probstei n  and 
Elliott,  and  the  Illingworth  transformation  sometimes  referred  to  as 
the  I  1 1 ingworth-Levy  transformation. 

In  the  Mangier  transformation  the  assumption  Is  made  that 
6  <<  rQ(x)  and  hence  r(x,y)  is  approximately  equal  to  rQ(x).  However, 
for  boundary  layer  thicknesses  of  the  order  of  r0(x),  this  approxi¬ 
mation  is  not  justified.  In  this  light  Probstein  and  Elliott  used  the 
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following  transformation  to  transform  the  axisymmetric  boundary  layer 
equations  containing  second-order  transverse  curvature  terms  to  an  almost 
two-dimensional  form  in  the  independent  variables  x  and  y 

dx  =  r.°2^L)..  dx  CIO) 

L2 


dy  =  r.(x>Y±  dy  (11) 

This  transformation  differs  from  Tangier’s  original  transformation  in 
that  r(x,y)  is  used  to  replace  rQ(x)  in  the  change  of  the  independent 
variable  y. 

In  working  with  the  boundary  layer  equations  in  two-dimensional 
form,  Illingworth  C6)  used  a  transformation  that  has  been  expressed  by 
Levy  (7)  in  the  following  form 


x  ___ 

M  =  /  PePeue  dx  < 12) 

o 


Combining  Equations  10  and  II  with  Equations  12  and  13,  redefining  the 
limits  of  integration,  and  introducing  two  new  variables  $  and  n  to 
replace  M  and  N,  yield  the  following  transformation  variables 


5  =  I  *Jepeue  X  dx 
o  L2 


■  9 


(14) 
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H  = 


(15) 


It  might  be  pointed  out  that  the  transformations  can  be  carried  out  with¬ 
out  introducing  the  characteristic  length  L. 


2.2.2  Transformed  Equations 

Applying  the  chain  rule  of  partial  differentiation  to  the 
transformation  variables,  the  necessary  operators  become 


»e  P 

*2 r 


r_ 

L 


Id  Peue^e  roZ 

[f>  1 

Jan 

3  l  ' 

My - P 

JjtJ 

n  N 

y 

KJ 

( 16) 


(17) 


Although  the  term  [lUi  does  not  appear  in  the  final  equations,  it  is 

is«iy 

included  here  for  completeness 


H 


Je 


'25 


}  i  3(Pr) 
J  L  35 


I  duc 


dy  +  dT 


-  f  _pJL_ 

Ul  J  L/2C 


dy 


(18) 


The  continuity  equation  is  identically  satisfied  by  intro¬ 
ducing  the  stream  function  <•> 


3  ip 

P™  =  3y 


(19) 
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prv  =  -  -i. 

3x 


(20) 


Applying  the  operator  given  by  Equation  16  to  Equation  19  gives 


'3* 

3n 


=  m  l  h- 


(21) 


.  1C  e 

Defining  a  non-dimensional  stream  function  f(£',n)  by 


*(C,n)  =  /2T  L  f(5,n) 


(22) 


one  has 


—  =  f'U,n>  (23) 

ue 

where  the  prime  denotes  differentiation  with  respect  to  n. 

Substituting  the  steady  flow  Euler's  equation 


dp  diie 

dx  pe  ue  — dx 


(24) 


into  the  momentum  equation  and  applying  the  operators  given  by  Equations 
16  and  17  to  Equations  2  and  5,  yield  the  following  forms  of  the 
momentum  and  energy  equations 


where 
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6  = 


.  2?  du* 


e 


(27) 


c  =  -£ H_  (28) 

P^e 


It  Is  of  Interest  here  to  consider  the  term  r2/rQ2  In  Equa¬ 
tions  25  and  26.  If  the  approximation  r  *  rQ  is  made,  as  for  thin 
boundary  layers,  the  term  r2/rQ2  is  one  and  the  resulting  equations  de¬ 
scribe  boundary  layers  without  second-order  transverse  curvature.  For 
this  case  the  transformation  given  by  Equation  15  must  contain  r0  rather 
than  r.  These  equations  can  also  be  used  to  describe  two-dimensional 
boundary  layers  by  setting  r2/rQ2  equal  to  one  In  Equations  25  and  26. 

For  this  case  the  transformations  given  by  Equations  14  and  15  must  have 
r02/ L2  and  r/L  set  equal  to  one.  Hence,  it  would  be  convenient  to 
Incorporate  Into  one  program  the  capability  of  handling  two-dimensional 
or  axisymmetric  boundary  layers  with  or  without  second-order  transverse 
curvature  terms.  This  is,  in  fact,  what  Jaffe,  Lind,  and  Smith  have  done. 

The  term  r2/r02  can  be  written  in  terms  of  the  variables  {■ 

and  n  as 


2L  /2T 
ro2Pe 


cos  a 


e 


Def i n i ng 


(30) 
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2L  /2g  cos  a  ^ 

ro2  Pe  ue  Jo 


(31) 


one  can  see  from  the  discussion  above  that  the  effects  of  transverse 
curvature  depend  upon  the  parameter  t.  Equations  25  and  26  now  become 


r 

» 

|C(l-t)fn 

+  8 

P 


+  ff' 


2£  f '  lH 
,  K 


(32) 


[0(|-+)  [f*F+ 


-  w 


t  fg»  =  2K  f' 


M  _  Q,  11 
H  9  H 


(33) 


with  the  boundary  conditions 
at  n  =  0 

f  =  0 

f'  =  0 

9  =  9w 

as  n  ■*  ne 

f*  -*■  I 

g  -*■  i 


(34a) 

(34b) 

(34c) 


(34d) 

(34e) 


2.2.3  Inverse  Transformations 

The  inverse  transformation  of  mapping  points  in  the  Cn-plane 
back  into  the  physical  xy-plane  usually  does  not  pose  a  difficult 
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problem,  especially  for  external  flow.  However,  for  internal  flow,  one 
should  determine  whether  the  inverse  transformation  exists  and  is  one- 
to-one  . 

An  existence  theorem  for  Inverse  transformations,  as  taken 
from  Olmsted  (8),  states  that  if:  (i)  £  and  n  are  single  valued  func¬ 
tions,  (ii)  £  and  n  are  continuously  differentiable  in  a  neighborhood  of 

the  point  (c,d),  and  (iii)  the  Jacobian  J  =  j-llzH)—  is  nonzero  at 

3(x,y) 

(c,d),  then  there  exists  a  neighborhood  about  the  image  of  the  point 
(c,d)  such  that  the  inverse  transformation  is  assured  and  the  corre¬ 
spondence  of  the  points  about  (c,d)  and  their  image  in  the  xy-plane  is 
one-to-one.  For  the  transformation  variables  given  by  Equations  14  and 
15,  the  only  difficulty  that  might  occur  is  due  to  the  Jacobian  J.  From 
the  definition  of  the  Jacobian  one  finds 

J  -  ..".e.'WP  r°2r  (35, 

lr 

There  are  two  places  where  the  Jacobian  is  zero.  One  is  In  the  plenum 
chamber  of  the  nozzle  where  ue  =  0.  This  point  will  be  discussed  in  the 
section  on  the  Solution  of  the  Transformed  Equations.  The  other  is 
along  the  centerline  of  the  nozzle  where  r  =  0.  This  point  could 
obviously  be  avoided  by  not  integrating  quite  all  the  way  to  the  center- 
line;  but,  in  the  numerical  solution  of  this  problem  no  difficulties 
were  encountered  in  obtaining  inverse  transformations  at  the  nozzle 
centerline.  Even  though  the  singularity  posed  no  problem  in  this  work. 
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the  point  to  be  made  is  that  in  working  with  internal  flows,  the  center- 
line  may  pose  certain  problems  that  require  special  attention.  This 
can  be  seen,  for  example,  in  the  work  of  Adams  (9). 

Only  the  inverse  transformation  of  n  is  necessary  here  since  E, 
is  calculated  for  each  prescribed  value  of  x.  The  inverse  of  Equation 
15  is  given  by 


y  = 


t m 

peue 


n 

r  fe  L 

J  o  r 


dn  +  n(c) 


(36) 


Applying  the  condition  that  y  =  0  at  n  =  0  gives  £1(£)  s  0  and  hence 


y  = 


(37) 


2.3  BOUNDARY  LAYER  PARAMETERS 

Certain  parameters  which  are  convenient  in  describing  the 
effects  produced  by  boundary  layers  are  the  displacement  thickness, 
momentum  thickness,  skin  friction,  heat  transfer,  and  Stanton  number. 

For  calculation  purposes  these  quantities  are  expressed  in  terms  of  the 
transformation  variables  Z,  and  n. 

The  displacement  thickness,  6*,  which  is  a  measure  of  mass- 
flow  defect,  is  defined  by 

6*  Ye 

|  2n  r  peue  dy  =  J  2irr  (pe  ue  -  pu)  dy  (38) 

o  o 
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which  takes  into  account  transverse  curvature.  The  upper  limit  of 
Integration  In  the  right  hand  side  of  this  equation  indicates  the  dis¬ 
tance  from  the  wall  that  the  numerical  Integration  has  been  carried  to 
satisfy,  within  prescribed  accuracy,  the  boundary  condition  u/u0  -*■  I. 

By  substituting  Equation  4  into  Equation  38,  the  left  hand  side  of  this 
equation  can  be  integrated  In  closed  form  yielding  a  quadratic  express¬ 
ion  for  6*.  The  two  solutions  for  6*  give  one  value  less  than  rQ  and 
one  greater  than  rQ.  Of  course,  the  one  applicable  to  internal  flow  is 
6*  less  than  rQ.  Using  Equation  37  to  express  the  result  in  the  trans¬ 
formed  plane  gives 


where  ne  Is  the  value  of  n  required  to  satisfy,  within  prescribed 
accuracy,  the  boundary  condition  u/u_  ■*  I. 

w 

The  momentum  thickness,  0,  which  is  a  measure  of  momentum-flow 
defect.  Is  defined  by 

/6  2nr  Peue2  dy  =  2vr  p(ueu  -  u^)  dy  (40) 

o  o 

Treating  this  similarly  to  6*  yields 

2L  /2£~  cos  a 
r02  pe  u 0 

The  skin  friction  coefficient  is  defined  by 
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lW 


h  Pe  ue/ 


(42) 


where  the  shear  stress  tw  is 


Tw  =  Pw  lf£ 


(43) 


w 


Applying  the  operator  given  by  Equation  16,  evaluated  at  the  wall,  to 
tw  gives  the  following  result  for 


2 

cf  =  f 


h  ro  Pw 


—  tl  f" 

L  p  Pw  Tw 


The  heat  transfer  at  the  wall  Is  given  by 


(44) 


-«w  •  k  If  J 


Treating  this  similarly  to  Cf  yields 


_  _  Pw  pw  ^o  ue  ro  9 1 w 

-  qw  - 


Prw  nr  i 

The  Stanton  number  is  defined  as 


(45) 


(46) 


St  = 


-c!w 


pe  ue  Hq  ( |-gw) 
which' in  the  transformed  plane  is 


(47) 


St  = 


V¥  Pw  ro  9'w 


(48) 


PrwL  pem  ( l-9w) 

Equation  37  can  be  used  to  obtain  y  In  the  form 


17 


AEDC-T  R-68-193 


Y  = 


cos  a 


I  -  (I  -  t)1* 


(49) 


where  +  Is  the  transverse  curvature  term  defined  by  Equation  31. 


2.4  SOLUTION  OF  THE  TRANSFORMER  EQUATIONS 


A  discussion  of  the  numerical  Integration  techniques  used  to 
solve  the  momentum  and  energy  equations  for  external  flow  can  be  found 
In  Reference  4.  Hence,  attention  here  will  be  focused  on  the  applica¬ 
tion  of  these  techniques  to  internal  flow. 

The  nozzle  wall  radius  rQ  Is  known  for  each  point  z  along  the 
axis  from  the  nozzle  geometry.  However,  the  coordinate  system  requires 
rQ  to  be  known  for  each  point  x  along  the  wall.  By  using  the  integral 
form  of  arc  length  one  can  write 


(50) 


If  the  nozzle  geometry  is  simple,  for  example  a  conical  nozzle  with  the 
converging  portion  described  by  an  arc  of  a  circle,  Equation  50  can  be 
integrated  in  closed  form  and  the  inverse  taken  to  give 


z  =  A(x) 


(51) 


Then  z  can  be  substituted  Into  the  known  expression 

rQ  =  B(z)  (52) 
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to  give 


rQ  =  B[A(x)]  =  rQ(x)  (53) 

If  the  nozzle  wall  cannot  be  expressed  accurately  by  simple  functions 
which  permit  integration  in  closed  form,  one  must  resort  to  numerical 
integration  and  tabulate  rQ  for  various  x  values. 

The  solutions  are  started  at  the  entrance  of  the  converging 
portion  of  the  nozzle  where  ue  is  zero.  The  total  conditions  in  the 
plenum  chamber  are  therefore  the  initial  conditions.  As  pointed  out 
previously,  one  has  a  problem  at  the  start  in  obtaining  y  since  ue  =  0. 
Recalling  the  expression  for  y  given  by  Equation  49,  one  notices  that  it 
contains  the  transverse  curvature  term  t.  The  wall  slope  is  undefined 
at  this  point,  and  t  is  zero.  Therefore,  Equation  49  gives  y  =  0.  This 
problem  is  alleviated  by  using  the  second  station  values  of  E  and  ue  in 
the  expression  for  t,  Equation  31. 

The  solution  at  the  first  station  is  obtained  by  assuming  a 
linear  temperature  distribution.  Also,  all  the  derivatives  in  the 
transformed  streamwise  direction  are  taken  as  zero  for  the  first  station. 

The  basic  inputs  to  the  program  are  the  total  pressure,  total 
enthalpy,  wall  temperature  distribution,  nozzle  geometry,  Prandtl  num¬ 
ber,  and  some  initial  estimated  pressure  distribution  along  the  axis. 
Since  a  given  nozzle  will  establish  its  own  pressure  distribution  for 
each  set  of  plenum  chamber  conditions,  this  pressure  distribution  must 
be  obtained  as  part  of  the  solution.  This  is  accomplished  by  iteration 
in  the  following  way.  From  the  displacement  thickness  calculated  at 
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each  station  an  effective  inviscid  area  is  calculated  using  the  corre¬ 
sponding  wall  radius.  Taking  the  ratio  of  this  area  to  the  throat  area 
a  new  pressure  is  calculated  at  each  station  from  one-dimensional  ex¬ 
pansion  theory.  The  resulting  pressure  distribution  is  then  used  as 
input  to  the  program  and  the  process  is  repeated  until  the  pressure  dis¬ 
tribution  converges  within  prescribed  accuracy.  As  shown  later,  the 
solutions  have  converged  after  two  or  three  iterations. 

It  should  be  pointed  out  that  the  throat  area  used  to  calculate 
the  pressure  distribution  upstream  of  the  throat  is  given  by  the  actual 
geometric  throat  area  for  the  first  iteration,  and  then  by  the  effective 
Inviscid  throat  area  calculated  from  the  previous  iteration  for  each 
successive  iteration.  All  pressure  distribution  calculations  downstream 
of  the  throat  use  the  effective  inviscid  throat  area  corresponding  to 
the  particular  iteration. 

All  the  solutions  calculated  were  for  nitrogen  since  it  was 
the  test  gas  used  in  the  particular  nozzles  considered.  In  order  to 
compare  some  of  the  solutions  with  experimental  data,  the  initial  con¬ 
ditions  used  correspond  to  actual  plenum  chamber  conditions. 

The  expression  taken  for  viscosity  in  all  the  calculations 
was  Sutherland's  Law,  and  the  Prandtl  number  used  was  Pr  =  0.7068. 
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SECTION  III 

RESULTS  OF  CALCULATIONS 

Various  solutions  were  obtained  for  four  different  nozzles 
which  are  in  operation  in  low  density  test  facilities  at  the  Arnold 
Engineering  Development  Center  (AEDC).  The  solutions  were  obtained  on 
a  CDC  1604  computer  in  the  von  KSrmSn  Facility  at  AEDC.  The  dimensions 
of  the  nozzles  are  given  in  Table  I. 


TABLE  I 

NOZZLE  DIMENSIONS 


Nozzle 

La 

mm 

Z 

rE 

R 

Type 

Mach  Three 

67.40b 

5.333 

64.00 

*  15.00 

4.425 

Con  i  ca  1 

Mach  Nine 

17.34 

.07343 

16.59 

*'2.407 

1.0 

Contoured 

Mach  Ten 

19.62 

.07405 

18.67 

3.833 

1.0 

Contoured 

Mach  Eighteen 

59.80 

.100 

57.40 

14.01 

1.185 

Con i ca 1 

aTerms  are  defined  in  the  nomenclature 
bAI  I  dimensions  are  in  Inches 


3.1  SOLUTIONS  FOR  MACH  THREE  NOZZLE 

The  Mach  three  nozzle  is  a  9.35°  half -angle  conical  nozzle 
which  operates  in  the  ARC  8V  Vacuum  Chamber  in  the  Aerospace  Environ¬ 
mental  Facility  at  AEDC.  The  nozzle  wall  Is  cooled  with  liquid  nitrogen 
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at  a  temperature  of  77#K,  to  reduce  the  boundary  layer  growth.  Due  to 
the  cold  wall,  this  nozzle  can  be  operated  at  various  plenum  chamber 
conditions  without  the  boundary  layer  merging.  The  plenum  chamber 
temperatures  range  from  290°K  to  I000°K  and  the  plenum  chamber  pressures 
range  from  100  to  1500  microns  of  mercury.  These  conditions  permit  a 
Mach  number  range  of  2.7  to  3.5  and  a  Reynolds  number  per  foot  range  of 
100  to  3600.  Solutions  were  obtained  for  plenum  chamber  conditions  of 
300°K  and  500  microns  of  mercury  .with  four  different  wall  temperature 
distributions. 

Figure  2  indicates  the  convergence  of  the  solutions  to  a 

certain  displacement  thickness.  Each  Iteration  represents  a  new  solution 

using  the  pressure  distribution  given  by  the  previous  iteration.  This 

solution  was  obtained  using  the  actual  nozzle  wall  temperature  dlstrl- 

« 

butlon  which  is  denoted  as  Tw  =  VI.  Wall  temperature  distributions  used 
for  the  Mach  three  nozzle  solutions  are  given  In  Figure  3. 

The  effect  of  wall  cooling  Is  shown  in  Figure  4  by  the  results 
of  the  constant  wall  temperature  solutions  of  I00°K  and  200°K.  The  dis¬ 
placement  thickness  is  more  Influenced  by  cooling  than  is  the  boundary 
layer  thickness.  This  significant  reduction  in  displacement  thickness  is 
due  to  the  higher  gas  density  near  the  wall. 

Based  on  the  indicated  effects  of  wall  cooling  in  Figure  4, 
solutions  were  obtained  to  Investigate  the  advantage  or  disadvantage  of 
cooling  the  nozzle  exit  region  more  thoroughly.  The  nozzle  is  presently 
cooled  upstream  of  the  throat  by  liquid  nitrogen  in  contact  with  the 
wall,  and  downstream  of  the  throat  by  liquid  nitrogen  pumped  through 
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Fig.  2 


Plot  of  Displacement  Thickness  for  Seecessive  Iterations 
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Fig.  4  Plot  of  Displacement  Thickness  and  Boundary  Layer  Thickness 
for  Two  Different  Wall  Temperatures  for  Mach  Three  Nozzle 
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copper  tubing- which  is  wrapped  around  the  nozzle  wall.  From  thermocouple 
measurements  taken  along  the  nozzle  wall  the  actual  temperature  distri¬ 
bution  is  given  by  Tw  =  VI.  It  is  anticipated  that  if  the  last  12  Inches 
of  the  nozzle  were  cooled  by  direct  contact  with  liquid  nitrogen,  the 
resulting  wall  temperature  distribution  would  be  as  shown  by  Tw  =  V2  of 
Figure  3,  page  24. 

The  Mach  number  distributions  along  the  nozzle  ax’is  for  the 
two  wall  temperature  distributions,  VI  and  V2,  are  given  in  Figure  5. 

The  Mach  number  at  the  exit  of  the  nozzle  is  Increased  by  the  Tw  =  V2 
distribution  from  3.29  to  3.40.  This  produces  a  decrease  in  Reynolds 
number  per  foot  from  1324  to  1225.  However,  the  axial  Mach  number 
gradient  is  increasedat  the  exit  plane  as  shown  in  Figure  5.  It  might 
be  that  neither  the  Tw  =  VI  nor  Tw  =  V2  case  provides  a  flow  suffi¬ 
ciently  free  of  axial  gradients  for  some  test  purposes.  In  this  light 
it  might  be  desirable  to  control  the  wall  temperature  in  order  to  pro¬ 
duce  a  displacement  thickness  distribution  that  will  eliminate  axial 
gradients  in  the  test  region. 

The  results  of  the  solutions  for  various  wall  temperature 
distributions  are  compared  in  Figures  6  through  13.  In  Figure  6  the 
velocity  distribution  for  Tw  =  V2  was  not  sufficiently  different  from 
that  for  Tw  =  I00°K  to  warrant  a  separate  curve.  The  boundary  layer 
thickness  appears  to  be  about  the  same  at  the  exit  whether  the  nozzle 
is  cooled  all  the  way,  or  only  in  the  region  of  the  converging  section 
and  the  exit.  Of  course  the  problem  of  axial  gradients  exists  as 
previously  pointed  out. 
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Fig.  6  Plot  of  Velocity  Distributions  at  Mach  Three  Nozzle  Exit 
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Fig.  7  Plot  of  Total  Enthalpy  Distribution: 
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Fig,  8  Piet  of  Static  Temperature 
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v/v£ 


Fig.  9  Plot  of  Reynolds  Number  Distributions  ot  Mach  Three  Nozzle  Exit 
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Fig.  10  Plot  of  Mach  Number  Distributions  at  Mach  Three  Nozzle  Exit 
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Fig.  12  Normalized  Momentum  Thickness  Distributions  in  Mach  Three  Nozzle 
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The  Reynolds  number  profile  given  in  Figure  9  indicates  that  a 
relatively  constant  Reynolds  number  region  may  exist  that  is  larger  than 
the  isentropic  core  region.  This  is  due  to  the  slight  "dip"  in  the  static 
temperature  profile  near  the  edge  of  the  boundary  layer  which  compen¬ 
sates  for  the  velocity  dropping  off.  Such  a  situation  may  be  desirable 
for  testing  purposes  when  the  usable  invlscid  flow  region  becomes  small. 

3.2  SOLUTIONS  FOR  MACH  NINE  NOZZLE 

The  Mach  nine  nozzle  is  a  contoured  water  cooled  nozzle  de¬ 
signed  by  the  method  of  Potter  and  Durand.  It  operates  In  Tunnel  L  of 
the  von  Kdrmdn  Facility.  This  nozzle  was  designed  for  a  Mach  number  of 
9.0  at  plenum  chamber  conditions  of  T0  =  2365°K  and  PQ  =  30.0  psia. 

From  the  work  of  Kins  low  and  Miller  (10),  It  has  been  con¬ 
cluded  that  the  vibrational  modes  freeze  somewhere  upstream  of  the 
throat  for  this  particular  nozzle  and  conditions.  Therefore,  effective 
values  of  T0  and  PQ  were  selected  to  yield  the  actual  gas  properties  in 
the  free-stream  assuming  frozen  flow  throughout.  The  values  used  were 
T0  =  2475°K  and  PQ  =  30.9  psia. 

From  the  assumption  of  free-stream  frozen  flow,  one  might  be 
concerned  about  the  possibility  of  vibrational  relaxation  existing  in 
the  portion  of  the  boundary  layer  near  the  wall.  Such  a  situation 
might  occur  due  to  lower  velocities  in  this  region  permitting  more 
molecular  collisions  per  unit  length  and,  thus,  relaxational  effects 
may  exist.  However,  it  will  be  shown  later  that  the  solution  given  by 
assuming  vibrational ly  frozen  flow  agrees  with  experiment  data. 
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Therefore,  it  is  concluded  that  if  vibrational  relaxation  exists  across 
the  boundary  layer,  It  is  negligible  for  this  particular  condition. 

Various  results  from  the  solution  for  the  Mach  nine  nozzle  with 
a  constant  wall  temperature  of  322°K  are  presented  In  Figures  14  through 
18.  This  wall  temperature  is  estimated  to  be  the  actual  one  existing  in 
the  nozzle.  An  approximately  constant  wall  temperature  occurs  due  to 
the  cooling  arrangement. 

The  effect  of  wall  temperature  on  Stanton  number,  momentum 
thickness,  displacement  thickness,  and  boundary  layer  thickness  is  shown 
in  Figures  19  and  20.  These  solutions  were  obtained  for  the  same  pres-, 
sure  distribution  and  two  different  constant  wall  temperatures  of  Tw  = 
I00°K  and  Tw  =  333°K.  The  displacement  thickness  is  reduced  about  the 
same  amount  as  the  boundary  layer  thickness.  However,  these  solutions 
are  not  indicative  of  the  actual  boundary  layer  that  would  exist  in  the 
nozzle  for  these  two  wall  temperatures  since  the  solutions  were  not 
Iterated. 

3.3  SOLUTIONS  FOR  MACH  TEN  NOZZLE 

The  Mach  ten  nozzle  is  also  a  water-cooled  contoured  nozzle 
designed  by  the  method  of  Potter  and  Durand  which  operates  in  Tunnel  L 
of  the  von  KdrmSn  Facility.  This  nozzle  was  designed  to  operate  at  a 
Mach  number  of  10,0  and  plenum  chamber  conditions  of  TQ  =  3090°K  and 
PQ  =  18.0  psia. 

The  work  of  Reference  10  indicates  that  the  flow  in  this 
nozzle  is  vibrational ly  frozen  downstream  of  the  throat.  The  effective 
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Fig.  14  Static  Temperature  Distribution  at  Mach  Nine  Nozzle  Exit,  TQ=  136°K 


38 


AEOC-TR-68-193 


Fig.  16  Mach  Number  and  Velocity  Distributions  at  Mach  Nine 
Nozzle  Exit,  H,/H0  =  0.13,  Me  =  9.263 
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Fig.  17  Skin  Friction  Coefficient  in  Mach  Nine  Nozzle,  Hw/H0  =  0.13 
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Fig.  18  Momentum,  Displacement,  and  Boundary  Layer  Thicknesses 
in  Mach  Nine  Nozzle,  Hw/H0  =  0.13 
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Fig.  20  Plot  of  Momentum,  Displacement,  and  Boundary  Layer  Thicknesses 
in  Mach  Nine  Nozzle  for  a  Prescribed  Pressure  Distribution 
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T0  and  PQ  used  with  the  assumption  of  frozen  flow  throughout  the  nozzle 
was  Tc  =  3I20°K  and  Pc  =  18.58  psla.  The  wall  temperature  was  taken  as 
constant  at  333°K. 

Solutions  for  this  nozzle  were  calculated  to  compare  with 
experimental  heat  transfer  measurements  made  by  Carden  (II).  Since  the 
heat  transfer  rate  varies  more  rapidly  near  the  throat,  a  smaller  step 
size  than  usual  was  taken  for  the  calculations  in  this  region.  Solutions 
were  iterated  to  about  two  inches  downstream  of  the  throat  where  the  heat 
transfer  rate  becomes  relatively  constant  and  rather  smali  compared  to 
the  throat  values.  The  step  size  taken  in  this  region  ranged  from  0.020 
inches  to  0.100  inches  depending  upon  the  anticipated  pressure  gradient 
at  a  particular  point.  It  was  then  increased  by  about  a  factor  of  two 
without  any  significant  change  in  the  results. 

The  boundary  layer  In  this  nozzle  grows  very  rapidly,  and 
therefore  provides  a  good  example  of  the  development  of  a  low  density 
boundary  layer  with  a  very  large  favorable  pressure  gradient.  Figures 
21  through  25  illustrate  the  development  of  the  boundary  layer  In  the 
throat  region.  Figure  22  indicates  that  the  Mach  number  at  the  throat 
is  nearly  constant  through  most  of  the  boundary  layer  due  to  the  nature 
of  the  velocity  and  temperature  distributions.  Notice  in  Figure  25  that 
negative  displacement  thicknesses  were  obtained  through  the  converging 
region  to  about  0.25  Inches  downstream  of  the  throat.  The  distributions 
in  Figure  26  are  the  results  of  the  initial  input  pressure  distribution 
since  successive  iterations  were  not  obtained  in  this  region.  However, 
it  should  be  indicative  of  the  type  of  distributions  which  would  exist 
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Fig.  21  Velocity,  Mach  Number,  and  Static  Temperature  Distributions  at 
0.16  Inches  Upstreom  of  Mach  Ten  Nozzle  Throat,  Me  =  0,459 
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Fig.  23.  Velocity,  Mach  Number,  and  Static  Temperature  Distributions  at  0.51 

Inches  Downstream  of  Mach  Ten  Nozzle  Throat,  M  =  3.331 
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Fig.  24  Velocity,  Mach  Number,  and  Static  Temperature  Distributions  at  1.57 

Inches  Downstream  of  Mach  Ten  Nozzle  Throat,  M  =  5.094 
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Fig.  25  Momentum,  Displacement,  and  Boundary  Layer  Thicknesses 
in  Throat  Region  of  Mach  Ten  Nozzle 
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Fig.  26  Plot  of  Indicated  Distributions  at  10.20  Inches 
Downstream  of  Mach  Ten  Nozzle  Throat 
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across  the  boundary  layer  farther  downstream. 

The  plane  at  1.57  inches  downstream  of  the  throat  is  located 
at  only  14  per  cent  of  the  total  axial  distance  of  the  nozzle,  but  a 
Mach  number  of  5.094  has  been  reached  which  is  about  half  of  the  exit 
value  of  10.15.  The  velocity  at  this  point  has  attained  93  per  cent  of 
the  exit  value.  Such  an  expansion,  and  associated  pressure  gradient, 
provides  an  excellent  test  for  the  theoretical  solutions. 

3.4  SOLUTIONS  FOR  MACH  EIGHTEEN  NOZZLE 

The  Mach  eighteen  nozzle  is  a  water-cooled  14°  half-angle 
conical  nozzle  which  operates  In  Tunnel  M  of  the  von  Kdrmgin  Facility. 

The  plenum  chamber  conditions  used  for  this  nozzle  were  TQ  =  6I60°K  and 
PQ  =  203  psia,  and  the  wal I  temperature  was  taken  as  constant  at  300°K. 

Solutions  for  this  nozzle  were  iterated  twice  in  an  attempt 
to  ''bracket"  the  final  result.  One  might  suspect  that  such  a  thing  can 
be  done  by  referring  to  Figure  2,  page  23.  in  Figure  2  the  Initial 
Input  pressure  distribution  was  more  favorabie  than  actually  existed 
in  the  nozzle,  and  the  first  two  iterations  bracketed  the  final  result. 
By  initially  inputting  a  less  favorable  pressure  distribution  than  was 
expected  to  exist  in  the  nozzle,  it  is  anticipated  that  the  first  two 
iterations  would  again  bracket  the  final  result.  The  reason  for  such  an 
oscillation  Is  that  a  less  favorable  pressure  gradient  would  produce  a 
thinner  boundary  layer  and  then  the  new  pressure  gradient  calculated 
would  be  larger  than  the  first  and  hence  the  next  iteration  would  yield 
a  thick  boundary  layer.  This  oscillation  seems  to  occur  in  each  case; 
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however,  as  has  already  been  pointed  out,  the  convergence  is  very  rapid 
after  two  or  three  iterations. 

The  second  iteration  was  not  carried  out  completely  as  indi¬ 
cated  in  Figures  27  and  28  by  the  dashed  lines.  The  reason  was  that  the 
machine  time  required  for  this  second  iteration  was  too  long.  The 
initial  input  pressure  distribution  consists  of  many  closely  spaced 
values  permitting  a  smooth  distribution  and  hence  good  streamwise  deriv¬ 
atives.  With  a  smooth  distribution  of  properties  along  the  centerline  a 
relatively  large  step  size  could  be  taken  for  the  first  iteration.  How¬ 
ever,  this  caused  numerical  problems  in  the  second  iteration  since  the 
new  pressure  distribution  could  only  be  calculated  from  the  few  stations 
used  in  the  first  iteration.  The  streamwise  derivatives  were  then  not  as 
accurate  due  to  the  few  points  availabie  for  calculating  the  derivatives. 
Rather  than  use  more  stations  in  the  second  iteration  or  curve  fit  the 
calculated  pressures  in  order  to  smooth  +he  s+reamwise  derivatives,  the 
second  iteration  was  extrapolated  In  Figures  27  and  28  by  the  nature  of 
the  first  iteration.  Some  results  of  These  iterations  are  given  in 
Figures  27,  28,  and  29. 

3.5  COMPARISON  OF  SOME  OF  THE  SOLUTIONS 

As  shown  in  Figure  4,  page  25,  for  the  Mach  three  nozzle,  the 
effect  of  wall  cooling  is  to  produce  a  larger  relative  reduction  in  dis¬ 
placement  thickness  than  in  boundary  layer  thickness.  This  produces  a 
somewhat  surprising  result  in  that  a  nozzie  designed  for  a  certain  Mach 
number  and  throat  size  will  have  a  smaller  area  of  uniform  core  if  the 
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Fig.  27  Heat  Transfer  Rate  in  Mach  Eighteen  Nozzle 
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Fig-  29  Velocity  Distributions  from  First  and  Second  Iterations  of  Mach 
Eighteen  Nozzle  Solution  at  19.42  Inches  Downstream  of  Throat 
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nozzle  is  designed  for  a  cold  wail  ra+her  ihan  a  hot  wall.  The  reason 
for  this  is  that  the  nozzle  wail  must  be  changed  for  each  design  rather 
than  the  effective  Invisc'd  radius,  and  Therefore,  the  size  of  the  uni¬ 
form  core  is  less  for  the  cold  wail  since  the  difference  between  the 
displacement  and  boundary  layer  thickness  is  greater.  Of  course,  the 
actual  nozzle  radius  for  the  hotter  walf  wi!l  be  iarger,  and  hence  the 
wall  angle  will  be  greater  since  the  throat  in  each  case  is  the  same 
size.  Therefore,  the  limitation  of  designing  nozzles  with  hot  walls  is 
governed  by  the  permissible  physical  size  ot  the  nozzle  and  The  maximum 
expansion  angle  the  flow  can  stand  without  separating.*  Potter  and 
Carden  give  an  example  of  the  design  of  a  given  Mach  number  nozzle  which 
clearly  illustrates  a  reduction  In  uniform  core  due  to  using  a  cold  wall. 

However,  if  one  is  Interested  in  increasing  the  size  of  the 
uniform  flow  of  a  given  nozzle,  extreme  wall  cooling  is' an  advantage.  As 
pointed  out  earlier,  a  significant  axiai  gradient  may  aiso  occur  which 
might  be  alleviated  somewhat  by  a  variable  wall  temperature. 

It  is  necessary  In  using  integral  techniques  to  assume  an 
appropriate  velocity  profile.  Some  of  the  more  simple  analytical  ex¬ 
pressions  used  for  the  velocity  are  given  and  plotted  in  Figure  30.  The 
exponent  x,  has  been  altered  slightly  In  tne  expression 


u 

u£ 


=  sin  yt 


(54) 


in  an  attempt  to  get  better  agreement  and  still  keep  the  expression 
simple.  The  boundary  layer  thickness  used  here  ;s  that  defined  as  the 
point  where  u/uQ  =  0,995.  In  Figure  31,  better  agreement  with  the 
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expression  given  by  Equation  54  was  obtained  by  defining  the  boundary 
layer  thickness  as  the  point  where  u/uQ  =  0.999.  Since  the  velocity  u 
approaches  uQ  asymptotically,  there  can  be  a  significant  difference 
between  the  results  of  the  two  definitions  of  boundary  layer  thickness. 

It  should  be  noted  that  if  6  is  defined  by  u/ue  =  0.999  for  the  Mach 
three  nozzle  case  given  In  Figure  30,  the  resulting  agreement  with  the 
Indicated  expressions  is  worse.  The  velocity  profiles  in  the  Mach  three 
nozzle  all  appear  to  have  a  near  linear  portion  near  the  wall,  whereas 
all  other  profiles  appear  to  have  slightly  more  curve  to  them  near  the 
wall. 

t 

The  streamwise  distributions  of  (f*')w,  (g?)w,  and  $  are  given 
in  Figures  32  through  34.  The  x/L  value  corresponding  to  a  particular 
point  in  the  nozzle,  such  as  the  throat,  depends  on  the  particular  nozzle 
considered.  However,  the  relatively  sharp  ’'peak”  of  each  curve  in 
Figures  32  through  34  occurs  at  the  throat  of  the  indicated  nozzle. 

Numerical  difficulties  were  sometimes  encountered  relatively 
far  downstream  near  the  edge  of  the  boundary  layer  when  the  boundary 
layer  approached  the  thickness  of  the  nozzle  radius.  This  is  believed 
to  be  due  to  the  transformation  variables  used.  Consider  Equation  15 
which  is  the  transformation  applied  to  the  independent  variable  y.  At  a 
given  x  location,  everything  in  this  equation  is  constant  except  p,  r, 
and  y.  Moreover,  everything  in  this  equation  except  p,  r,  and  y,  is 
relatively  constant  for  all  x  values  near  the  exit  of  the  nozzle  since 
uQ  and  5  change  very  little  in  this  region.  However,  p  continues  to 
decrease,  and  therefore,  a  larger  upper  limit  of  integration,  y,  is 
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Fig.  33  Distribution  of 
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Fig.  34  Distribution  of  /3  in  Eoch  of  the  Four  Nozzles  Considered 
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required  in.  Equation  i5  to  satisfy  a  prescribed  value  of  n.  Approxi¬ 
mating  Equation  15  by  writing  it  in  finite  difference  form,  one  can  see 
that  a  given  increment  of  n,  which  Is  what  the  program  uses  in  the 
numerical  integration  across  the  boundary  layer,  requires  an  increasing 
increment  in  y  as  the  centerline  is  approached.  This  increasing  incre¬ 
ment  in  y  Is  due  to  the  fact  that  the  variable  r  is  small  near  the  center¬ 
line.  This  increment  of  y  is  sometimes  too  large  and  convergence  is 
difficult  to  obtain.  Hence,  the  centerline  Introduces  a  problem  which 
does  not  exist  in  external  boundary  layer  analysis  using  the  same 
transformation  variables. 
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SECTION  IV 

COMPARISON  OF  THEORETICAL  RESULTS  AND  EXPERIMENTAL  DATA 

It  Is  not  possible  to  compare  the  solution  of  all  the  dependent 
variables  with  corresponding  data  since  appropriate  experimental  data 
cannot  be  obtained.  However,  the  pitot  pressure  can  be  obtained  experi¬ 
mentally  and  it  provides  a  somewhat  satisfactory  test  of  the  theory  since 
It  is  dependent  upon  both  the  velocity  and  temperature.  Carden  has  made 
some  nozzle  wall  heat  transfer  measurements  which  can  be  applied  to  the 
throat  region  of  the  Mach  ten  nozzle.  Also,  Potter  and  Carden  present 
both  pitot  pressure  measurements  and  relative  heat  flux  measurements 
taken  from  the  Mach  nine  nozzle. 

The  experimental  data  which  are  compared  to  theory  consist  of 
measurements  of  pitot  pressure  at  the  exit  of  the  Mach  three  nozzle, 
measurements  of  pitot  pressure  and  relative  heat  flux  for  the  last 
three  inches  of  the  Mach  nine  nozzle,  and  measurements  of  heat  transfer 
rates  of  the  throat  region  of  the  Mach  ten  nozzle. 

Comparison  of  the  measured  pitot  pressures  with  the  calcu¬ 
lated  pitot  pressure  distribution  Is  given  In  Figure  35  for  the  Mach 
three  nozzle.  The  error  band  represents  the  accuracy  in  which  the 
actual  location  of  the  pitot  pressure  probe  was  known.  The  reason  for 
the  inaccuracy  Is  that  the  probe  was  mounted  on  a  remotely  control  led 
moveable  support  whose  position  was  only  known  within  t0.5  inches. 
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y/y^ 

Fig.  35  Plot  of  Calculated  and  Measured  Pitot  Pressure  Distributions 
at  Exit  of  Mach  Three  Nozzle  for  Tw  =  VI 
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These  pitot  pressure  measurements  were  corrected  for  viscous 
effects  by  using  the  calculated  Mach  number  and  Reynolds  number  corre¬ 
sponding  to  the  point  of  the  measurement.  The  correction  required  on  the 
four  points  farthermost  from  the  wall  ranged  from  one  to  five  per  cent. 
The  one  at  y/y^  =  0.07  required  43  per  cent  correction  and  the  one 
nearest  the  wall  was  nearly  a  static  pressure  measurement  and  no 
correction  was  applied. 

A  significant  result  of  the  comparison  in  Figure  35  is  that 
the  boundary  layer  equations  appear  to  adequately  describe  the  viscous 
region  in  a  highly  rarefied  flow.  For  instance,  the  mean  free  path  of 
the  gas  for  these  plenum  chamber  conditions  is  on  the  order  of  one-tenth 
of  an  inch  for  about  the  last  half  of  the  nozzle  length.  If  a  Knudsen 
layer  Is  defined  as  ten  mean  free  paths  from  the  nozzle  wall,  the  outer 
edge  of  the  Knudsen  layer  would  correspond  to  about  y/y^  *0.1.  This 
distance  is  about  25  per  cent  of  the  total  boundary  layer  thickness. 

Since  velocity  slip  and  temperature  jump  have  not  been  taken  into 
account,  the  agreement  is  surprisingly  good. 

A  comparison  of  calculated  and  measured  pitot  pressures  is 
given  in  Figure  36  for  the  Mach  nine  nozzle.  The  continuous  measured 
pitot  pressure  distribution  was  obtained  by  connecting  the  pitot  probe 
to  a  linear  potentiometer  which  in  turn  was  connected  to  an  x-y  plotter 
along  with  the  signal  from  the  pitot  probe  pressure  transducer.  Thus, 
it  was  possible  to  get  a  continuous  trace  for  the  pitot  pressure  as  a 
function  of  position  rather  than  discrete  points.  Some  experimental 
pitot  pressure  data  are  given  in  Figure  37  which  were  taken  from 
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Fig.  36  Plot  of  Calculated  and  Measured  Pitot  Pressure 
Distributions  at  Exit  of  Mach  Nine  Nozzle 
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Fig.  37  Plot  of  Calculated  Total  Enthalpy  Distribution,  Calculated  Edge  of  Boundary  Layer, 
Measurements  of  Relative  Heat  Flux,  and  Measurements  of  Pitot  Pressure  for  Mach 
Nine  Nozzle 
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Reference  2.  The  previously  mentioned  method  of  taking  pitot  pressures 
was  not  available  for  use  at  the  time  these  data  were  obtained. 

The  pitot  pressure  measurements  presented  In  Figures  36  and 
37  were  not  corrected  for  viscous  effects  since  the  correction  required 
for  the  measurements  over  most  of  the  boundary  layer  ranged  from  zero  to 
two  per  cent.  The  region  of  the  boundary  layer  near  the  wall  in  Figure 
36  probably  requires  a  significant  correction,  but  good  experimental 
data  for  viscous  effects  are  not  available  for  the  flow  conditions  in 
this  region,  so  a  correction  was  not  attempted.  However,  there  are 
quantitative  data  available  which  indicate  that  any  correction  applied 
to  the  measured  pitot  pressures  near  the  wall  would  lower  the  measured 
values  and  thus  yield  better  agreement  with  the  calculated  distribution. 
Therefore,  it  must  be  concluded  that  this  comparison  between  measured  and 
calculated  pitot  pressures  is  rather  good. 

A  pitot  pressure  survey  such  as  that  given  in  Figure  37  is 
usually  the  only  information  available  for  determining  the  boundary 
layer  thickness.  It  is  usually  a  good  assumption  to  consider  the  pitot 
pressure  distribution  approximately  the  same  as  the  velocity  distribu¬ 
tion  near  the  edge  of  the  boundary  layer.  There  are  situations,  however, 
where  this  assumption  is  not  justified.  Examples  of  this  are  illustrated 
in  Figures  21  through  24,  pages  46  through  49.  The  pitot  pressure  Is  a 
function  only  of  the  Mach  number  at  each  axial  location  since  the  static 
pressure  Is  assumed  constant;  therefore,  the  pitot  pressure  will  drop 
off  at  the  same  point  as  the  Mach  number.  For  these  cases  one  can  see 
that  a  pitot  survey  would  be  very  misleading  as  to  the  location  of  the 
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edge  of  the  boundary  layer  since  the  Mach  number  is  constant  through 
much  of  the  boundary  layer.  Fortunately,  near  the  exit  of  supersonic 
and  hypersonic  nozzles  where  most  of  the  total  energy  of  the  flow  has 
been  converted  into  velocity,  the  temperature  distribution  across  the 
boundary  layer  is  such  that  the  Mach  number,  and  hence  the  pitot  pres¬ 
sure,  is  usually  a  good  way  to  find  the  edge  of  the  boundary  layer. 
However,  for  subsonic  flows,  and  supersonic  flows  where  the  static 
temperature  is  high  relative  to  the  wall  temperature,  the  pitot  pres¬ 
sure  is  not  a  good  means  for  such  an  investigation. 

The  measured  pitot  pressures  in  Figure  37  were  used  to  estimate 
6  for  the  Mach  nine  nozzle.  The  calculated  values  of  6  are  indicated  for 
comparison  with  the  measured  pressures.  For  the  conditions  at  the  exit 
of  this  nozzle,  the  measured  values  of  6  are  consistent  with  the  cal¬ 
culated  values.  Table  II  gives  a  tabulated  comparison  of  the  calculated 
and  experimental  values  of  Mach  number  and  displacement  thickness  at 
the  exit  of  the  Mach  three  and  nine  nozzles.  The  actual  displacement 
thickness  was  calculated  using  the  acTual  Mach  number  and  actual  nozzle 
wall  radius. 


TABLE  II 

COMPARISON  OF  THEORY  AND  ACTUAL  VALUES 


"1 

w  -  '■ 

5*( i nches) 

Nozzle 

mi 

Calculated 

1 H  3531 ;  I 

Ca 1 cu 1 ated 

Mach  Three 

3.28 

3.29 

2.47 

2.41 

Mach  Nine 

9.30 

9.27 

0.95 

0.96 
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Again  there  is  consistent  agreement  between  the  experimental 
data  and  numerical  solutions.  In  Reference  li,  Carden  has  compared 
similar  solutions  with  these  experimental  heat  transfer  rate  measure¬ 
ments.  The  agreement  is  good  especially  downstream  of  the  throat.  How¬ 
ever,  the  non-similar  solutions  presented  herein  predict  the  heat  trans¬ 
fer  rate  at  the  throat  more  accurately  than  do  the  similar  solutions. 
However,  this  would  be  expected  since  the  conditions  for  similar 
solutions  are  more  difficult  to  satisfy  near  the  throat. 
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Fig.  38  Plot  of  Calculated  and  Measured  Heat  Transfer  Coefficients  in  Mach  Ten  Nozzle 
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SECTION  V 
CONCLUSIONS 

The  ax i symmetric  laminar  boundary  layer  equations  including 
second-order  transverse  curvature  terms  have  been  adapted  to  the  pur¬ 
pose  of  describing  internal  laminar  boundary  layers  associated  with 
axlsymmetric  nozzles.  These  equations  have  been  formulated  for  solu¬ 
tion  similar  to  the  method  used  by  Jaffe,  Lind,  and  Smith  for  the 
external  boundary  layer  equations.  This  approach  was  taken  in  order 
to  take  maximum  advantage  of  the  numerical  integration  techniques  which 
Jaffe,  Lind,  and  Smith  have  developed. 

Non-slmI!ar  solutions  were  obtained  for  four  existing  low 
density  nozzles  using  various  wall  temperature  distributions.  Appro¬ 
priate  solutions  were  compared  to  experimental  data  which  were  taken  for 
the  specific  purpose  of  testing  the  theoretical  results.  Also  com¬ 
parisons  were  made  with  data  that  have  been  taken  In  the  past  by  Potter 
and  Carden.  In  all  cases  the  theoretical  solutions  appeared  to  be  in 
substantial  agreement  with  the  experimental  data.  Therefore,  it  5s  con¬ 
cluded  that  this  set  of  equations,  which  Is  an  approximation  +o  the 
general  form  of  the  continuity  equation,  Navier-Stokes  equations,  and 
energy  equation,  adequately  describes  the  viscous  flow  region  in  the  low 

I 

density  axlsymmetric  nozzles  considered. 

Probably  the  most  important  result  of  the  investigation  is  the 
capability  to  analyze  the  boundary  layer  In  an  existing  or  hypothetical 


74 


AEDC-T  R -68-193 


ax i symmetric  nozzle  of  arbitrary  geometry,  plenum  chamber  conditions,  and 
wait  temperature  distribution.  The  satisfaction  of  such  a  need  Is  both 
economical  and  practical. 

The  boundary  layer  equations  have  been  shown  to  adequately 
predict  physical  phenomena  in  a  gas  that  is  highly  rarefied.  From  this 
result  it  would  seem  appropriate  to  apply  the  boundary  layer  equations 
to  the  slip  flow  and  temperature  jump  regime  by  properly  modifying  the 
wall  boundary  conditions. 

interesting  and  useful  extensions  of  this  work  would  be  the 
following:  (I)  modify  the  transformation  variables  and/or  the  program 
to  overcome  the  aforementioned  difficulty  near  the  nozzle  centerline, 

(2)  Include  velocity  slip  and  temperature  jump  In  the  wall  boundary 
conditions,  (3)  Include  the  possibility  of  non-equl I Ibrlum,  and  (4)  apply 
this  non-similar  method  to  the  design  of  nozzles.  Actually  this  latter 
extension  would  be  rather  simple,  since  the  convergence  of  the  solution 
could  be  accomplished  by  fixing  the  effective  invlscid  nozzle  wall  and 
iterating  on  the  actual  nozzle  wall  rather  than  vice  versa  as  it  is  now 
done,  in  this  case  the  Input  to  the  program  for  each  successive 
iteration  would  be  the  actual  nozzle  wall  radius  calculated  by  the 
previous  iteration.  The  properties  along  the  centerline  would  of  course 
be  inose  for  which  the  nozzle  is  Intended  to  achieve.  These  properties 
are  usually  obtained  from  an  Invlscid  method  of  characteristics  solution 
by  specifying  the  desired  plenum  chamber  conditions,  throat  radius, 
maximum  expansion  angle,  and  uniform  flow  conditions  at  the  nozzle  exit. 
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